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Abstract 

We determine tlie phase diagrams of conservative diffusive contact processes by means of numer- 
ical simulations. These models are versions of the ordinary diffusive single-creation, pair-creation 
and triplet-creation contact processes in which the particle number is conserved. The transition 
between the frozen and active states was determined by studying the system in the subcritical 
regime and the nature of the transition, whether continuous or first order, was determined by 
looking at the fractal dimension of the critical cluster. For the single-creation model the transi- 
tion remains continuous for any diffusion rate. For pair- and triplet-creation models, however, the 
transition becomes first order for high enough diffusion rate. Our results indicate that in the limit 
of infinite diffusion rate the jump in density equals 2/3 for the pair-creation model and 5/6 for the 
triplet-creation model. 

FAGS numbers: 05.50.Ln, 05.50+q, 05.65.+b 



1 



I. INTRODUCTION 

The usual contact process is a model for nonequilibrium process composed by 

two subprocesses: a catalytic creation and a spontaneous annihilation. In the usual contact 
process, which we call single-creation contact process, aparticle is created if the site has at 
least one neighbor site occupied. Many generalizations can be made by taking into 

account the basic mechanisms contained in the single-creation model. In the pair-creation 
contact process, for instance, at least two neighbor sites occupied are necessary to create 
a new particle. In the triplet-creation model one should have at least three sites occupied. 
All these variants of the contact process exhibit a continuous phase transition between an 
absorbing state and an active state that belongs to the direct percolation (DP) universality 
class. 

Diffusive models are defined by permitting the diffusion of particles in addition to the 
catalytic creation and spontaneous annihilation. A diffusion process is done by moving a 
particle to an empty nearest neighbor site. The introduction of diffusion does not destroy 
the transition from an active state to an absorbirig state giving rise to a transition line that 
separates the two phases. Jensen and Dickman [7| have shown that for the single-creation 
diffusive contact process this line is always second order for any diffusion rate and belongs to 
the DP universality class, i. e., the addition of the diffusion does not change the universality 
class nor the nature of transition. For the pair-creation the numerical results by Dickman 
and Tome show that the transition is continuous for low diffusion and they argue that 
there is no change in the nature of the transition for high diffusion. For the triplet-creation 
contact process, Dickman and Tome [s] have shown that the transition becomes a first order 
for high enough diffusion. 

In the present work, we study conservative versions of the models mentioned above. 
A conservative version of a model for nonequilibrium process was introduced by Ziff and 
Brozilow 91], who used a const ant- density ensemble to study the ZGB model. A conservative 
version of directed percolation was used by Broker and Grassberger |lO]. The conservative 
contact process was introduced by Tome and de Oliveira ll|] who have shown its equivalence 
in the thermodynamic limit to the ordinary contact process and how to calculate the rates 
from averages over the const ant- density ensemble. The equivalence between this ensemble 
and the constant-rate ensemble was later proved by Hilhorst and Wijland 121]. 
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In the conservative versions of contact processes ll|, [is], Q], an empty site becomes 
occupied in a way similar to the catalytic creation. But instead of creating a new particle, 
like in the ordinary contact processes, a randomly chosen particle of the system leaves its 
place and jumps into the empty site. Thus, both the creation and annihilation processes 
are replaced with a jumping process. However, this is not an unrestricted jumping because 
particles are not allowed to jump to a vacant site surrounded by empty sites. It is necessary 
to have a neighborhood (a set of one, two, or three sites, depending on the model) of sites 
occupied. 

One advantage of using the conservative versions is that they allow us to study the 
model without the danger of falling down into the absorbing state. The conservation of 
particles permits us to perform numerical simulations that avoid the accidental fall into the 
absorbing state. Although they do not have absorbing states, they are equivalent, in the 
thermodynamic limit, to the ordinary models. The conservative and ordinary models are 
versions of the same model in distinct ensembles 
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; the first models belonging to 



the constant-particle ensemble, the second models belonging to the constant-rate ensemble. 

Another advantage is related to the expected existence of first order transition. In the 
ordinary models, a very small change in the annihilation rate (the control parameter), near 
the transition, causes a great change in the density. In the conservative models, because 
of the fact that the particle number is a conserved quantity (and therefore, works as the 
control parameter), this problem does not occur. This advantage has been used by Ziff and 
Brosilow in their study first order transition in the ZGB model. 



II. CONSERVED REACTION DIFFUSION MODELS 

In the construction of conserved models we have to be concerned only with the reaction 
process since the diffusion process already conserves the particle number. The necessary con- 
dition to set up an equivalent conserved version of an arbitrary ordinary reaction-diffusion 
process in a lattice is that the reaction process be a sum of a creation subprocess and an 
annihilation subprocess [l^. This is always possible do realize because these two subpro- 
cesses are mutually excludent. If a site of the lattice is empty only creation is possible; if 
it is occupied, only annihilation is possible. Therefore, the transition rate Wi related the 
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creation-annihilation of a particle at site i can always be written in the form 

Wi = kcujf + kAUjf, (1) 

where the first term is related to the creation of a particle at site i and the second to 
annihilation of a particle at site i. The quantity cup vanishes if there is already a particle at 
site i and uf vanishes if site i is empty. The quantities kc and A;^ are the actual parameters 
of the ordinary model which we call amplitudes of the creation and annihilation rates, 
respectively. 

The conserved version is set up by replacing both the creation and annihilation subpro- 
cesses by a particle jump process i — » j with rate Wij = ufuf / L where L is the number 



of sites of the lattice. One can prove [ij] that a two site process defined by this transition 
rate is equivalent in the thermodynamic limit to the ordinary process. To see how this come 
about let us look at the total rate Y.i'^ij iii which particles jump to site j. In the ther- 
modynamic limit, the sum Y.i^t/L approaches, by the law of large numbers, the average 
{ujf) so that Wij = {ujf)uj'^ . By an analogous argument the total rate in which particles 
leave the site i is 'Wij = {uj^)ujf. The averages {uf) and (cuf ) act then as the amplitudes 
of creation and annihilation rates, respectively, what allows us to write down the following 
relation [l^ 

(2) 



kA ^ (4^ 
kc {u^f) ' 



between the amplitude rates of the constant-rate ensemble and averages determined in the 
constant-particle ensemble. 

For the model we study here particles are spontaneous annihilated do that is 1 if site 
i is occupied and if it is empty. Therefore J2i — where n is the number of particles. 
The quantity uj^ is if site j is occupied. Since creation is catalytic, this quantity depends 
also on the neighborhood of site j. For the single-creation model it equals half the number 
of nearest neighbor occupied sites. For the pair-creation model it equals half the number 
of pairs of nearest neighbor occupied sites. For the triplet-creation model it equals half the 
number of triplets of nearest neighbor occupied sites. It is convenient do define an active 
site as the site for which cup is nonzero. The number of active sites is defined by 



E<- (3) 
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We also define a quantity a as being the right-hand side of equation ([2]) so that, for the 
models studied here, 

a = (4) 

n 

where the averages are taken in the constant-particle ensemble. Usually one defines the 
ordinary reaction process so that the rate amplitudes are kc = 1 and kA = k. Therefore, 
according to relation ([2]), a coincides with the parameter k of the ordinary model as long as 
the average density of particles of the ordinary models equals the density of particle n/L in 
the conserved models. 

The rules of the reaction-diffusion processes we used are such that the diffusion occurs 
with probability D and the jump process with probability 1 — D. The quantity D and the 
diffusion rate D are related by 

D = ^. (5) 

1 + D ^ ^ 

III. EXACT AND MEAN-FIELD RESULTS 

The average number of active sites per site of the lattice equals the probabilities P(10), 
P(llO), and P(lllO) for the single-creation, pair-creation and triplet-creation models, re- 
spectively. Since the number o particles per site is the probability P(l) it follows that a is 
given, respectively, by a = P(10)/P(l), a = P(110)/P(l), and a = P(1110)/P(l), for the 
three models. In the limit of infinite diffusion rate the particles will be uncorrelated so that 
P(10) = P(1)P(0), etc. Taking into account that P(l) = p and P(0) = 1 - p, we get the 
following exact results for the active state, valid for D = 1, 



for the single-creation model. 



for the pair-creation model, and 



a = 1 - p, (6) 



a = p(l-p), (7) 



a = p2(l-p), (8) 



for the triplet-creation model. 

These results give a continuous transition for the single-creation model. For the pair- 
creation and triplet-creation models, on the other hand, they give a discontinuous transition 
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since p does not vanishes continuously as one increases a. The quantity a has a maximum at 
a nonzero value po of the density which is po = 1/2 for the pair-creation model and po = 2/3 
for the triplet-creation. The corresponding values of a are = 1/4 and = 4/27, 
respectively. Since there is no free energy from which we could decide at what point the 
jump in the density occurs one is tempted to use the maximum value of a (spinodal point). 
However, as we will see, our numerical results do not support this point of view. According 
to the numerical results, the discontinuity occurs at a smaller value of a. 

For the single creation the diffusion does not change the nature of the transition. Even 
at infinite diffusion rate the transition is continuous as the exact result (Ej) shows. The 
critical line on the diagram D versus a can be obtained by a mean-field approximation. By 
using a two-site mean-field approximation we get a relation between p and a which shows a 
continuous transition for all values of D and which recovers the exact result ([6]) when D = 1. 
The critical line obtained from this approximation is given by 

^ 2a — 1 , . 

showing that a — 1 as D ^ 1 in accordance with the exact result ([6]) and, as we will see, 
with numerical simulations. 

The exact results and ([8]) for the pair-creation and triplet-creation models cannot 
be used to infer that the transition will remain discontinuous for finite diffusion rate even 
if the rate is large. An indication that the transition is continuous at low diffusion and 
discontinuous for sufficiently large diffusion, giving rise to a tricritical point, comes from 
mean-field approximations which can be done at several levels |l5|. At the level of three 
sites the mean-field approximation locates the tricritical point of the pair-creation model at 
Dt = 0.032 and at = 0.1687. For the triplet-creation model it is necessary to use a higher 
order of approximation. At the level of four sites the tricritical point occurs at Dt = 0.017 
Sj. Although both results are in qualitative agreement with our numerical simulations they 
are very low when compared with the figures coming from the numerical simulations. 

IV. NUMERICAL SIMULATIONS 

We have simulated the conservative diffusive contact process in a one- dimensional lattice. 
The actual simulation is performed as follows. At each time step a particle is selected at 
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FIG. 1: The effective number of active sites per particle a as function of particle density p for 
single-creation conservative contact process for some values of probability D. 

random, say a particle at site and one of its neighboring sites is chosen randomly, say site 
j. If this neighboring site is empty then we decide which process to perform: the diffusion 
of particles, occurring with probability D' , or the creation-annihilation process, occurring 
with probability 1 — D' . If the diffusion process is chosen then the particle at i hops to the 
neighboring site j. If the creation- annihilation process is chosen then any another particle 
of the system, including the one at site is chosen randomly and placed at site j. In the 
case of the pair-creation or triplet-creation models, however, this only happens if the chosen 
particle i has at least one or two nearest neighbor occupied sites, respectively. The relation 
between the probability D' we use in the simulation and the actual probability of diffusion 



D is D' = 2D/{1 + D) [16|. This is so because we are choosing a particle from a list of 
occupied sites and then choosing with equal probability one of the neighboring sites to place 
the particle. 



A. Supercritical regime 

The simulation was performed using lattices with L sites and periodic boundary condi- 
tions. The averages of the appropriates quantities were obtained from a number of Monte 
Carlo steps ranging from 10^ to 3 x 10^, after discarding a sufficient number of steps to reach 
the stationary state. In Figs. [H [2]and[3l we show the particle density p = ra/L as a function 
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FIG. 2: The effective number of active sites per particle a as function of the particle density p for 
pair-creation conservative contact process for some values of probability D. The horizontal straight 
line at a = 0.222 was obtained by extrapolation. 

of a, calculated by using formula (jlj), for several values of the hopping probability D. We 
have used L = 10^ and varied the number of particles n. As expected, for high diffusion 
rate the curves approach the exact behaviors given by equations IQ, ([7]) and ([8]). 

For the diffusive single-creation contact process, the transition is found to be continuous 
for all values of D. Increasing the diffusion probability D the critical value of a increases 
towards the value 1 when D 1 as expected. For the diffusive pair-creation and triplet- 
creation models the phase transition is continuous for low diffusion becoming discontinuous 
for high enough diffusion. The tricritical point occurs at Dt = 0.965 ± 0.010 for the pair- 
creation contact process and Dt = 0.945 ± 0.005 for the triplet-creation as we shall see 
shortly. Figs. [2] and [3] seems to show that this is indeed the case. 

To compare the behaviors corresponding to the second and first order transitions, we 
simulated the pair-creation model at D = 0.5 and D = 0.995 for various values of the 
system size L ranging from 50 to 10^. For the former case, D = 0.5, the plot of p versus 
a, shown in Fig. |H shows a continuous transition. That the transition is continuous is 
confirmed by the the data collapse of the data shown in the inset of Fig. |H For D = 0.995, 
the plot of p versus a, shown in Fig. [5l displays a jump when L ^ oo increases. Results 
similar to those of Figs. H] and [5] are also found for the triplet-creation model. 
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FIG. 3: The effective number of active sites per particle a as function of the particle density p 
for triplet-creation conservative contact process for some values of probability D. The horizontal 
straight line at a = 0.115 was obtained by extrapolation. 
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FIG. 4: The effective number of active sites per particle a as function of the particle density p 
for several values of L in the supercritical regime for D = 0.5 for the pair-creation model. The 
inset show the scaling plot of y = L^^'^\ac — ot\ versus x = L^^'^p, using the DP critical exponents 
P = 0.277 and z/ = 1.097. 
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FIG. 5: The effective number of active sites per particle a as function of the particle density p for 
several values of L in the supercritical regime for D = 0.995 for the pair-creation model. 

B. Subcritical regime 



To simulate the system in the subcritical regime we consider an infinite lattice with a finite 
number n of particles. In practice we use a finite lattice and check whether a particle reaches 
the border. If a particle is about to reach the border we increase the size of the lattice. For 
a fixed value of D we have simulated the system for several values of ra, computing a by 
using (jlj). For each value of Z), the critical value etc was obtained in the limit n ^ oo by 
a linear extrapolation in 1/n. Using these results we have built the phase diagram in the 
plane D versus a, as shown in Figs. [6l [7] and [H The numerical values we have obtained 
for the transition line agrees very well with the results obtained previously for the ordinary 



models 



7 



When D — s> 1 the critical value of a approaches a limiting value cto- For the single- 
creation model ao = 1 as expected from the exact result ([6]) and also from the mean-field 
result ([9]). Assuming that the behavior of D around a = 1 is given by 



(10) 



we have found from the plot shown in the inset of Fig. [6] that = 4.03(3). Notice that the 
mean-field behavior, given by ([9]), predicts the value = 1. 

An important feature of the models studied here is the emergence of a fractal structure at 
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FIG. 6: Phase diagram for the single-creation conservative contact process. The star corresponds 
the value of ac = 1 in the limit D =1. The inset corresponds the Log-log of equation (jlOp . The 
transition from active to frozen state is always second-order. 
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FIG. 7: Phase diagram for the pair-creation conservative contact process. The star corresponds 
the value of ao = 0.222 in the limit D =1. The tricritical point (full circle) is located at at = 0.199 
and A = 0.965. 

the transition point, characterized by its fractal dimension. We have calculated the fractal 
dimension at the transition for each value of D. To this end we have simulated a system 
with n particles and determined the average distance R between the two particles located 
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FIG. 8: Phase diagram for the triplet-creation conservative contact process. The star corresponds 
the value of ao = 0.115 in the limit D = 1 The tricritical point (full circle) is located at at = 0.102 
and Dt = 0.945. 

at the extremities of the system. We assume the asymptotic behavior [l^ 

n ~ R^^, (11) 

where dp is the fractal dimension, so that the slope of a log-log plot of versus R gives 
the fractal dimension as shown in Figs. [9] and [TOl for the pair and triplet- creation models, 
respectively. 

For a continuous phase transition studied here we expect the emergence of a fractal cluster 
with a fractal dimension dp strictly less than one. Indeed for the single-creation model and 
for the pair and triplet up to a certain value of D we found a fractal dimension dp = 0.75 
which is the expected value for a system in the DP universality class. However, for the pair 
and triplet creation and for sufficient large values of the fractal dimension becomes the 
Euclidean dimension dp = d = 1 reflecting the formation of a compact cluster whose size 
R increases linearly with n. The changing of behavior occurs at Dt = 0.965 ± 0.010 and 
at = 0.199 ±0.003 for the pair-creation model and A = 0.945 ±0.005 and at = 0.102 ±0.001 
for the triplet-creation model. The tricritical point obtained by Dickman and Tome [8] for 
the ordinary triplet-creation model by means of numerical simulations is Dt ~ 0.85 and 
at — 0.096. These values correspond, actually, to a point over the critical line in Fig. [81 
For D = 0.85 our result is a = 0.0954. As for the pair-creation model, the numerical results 
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FIG. 9: Log-log plot of size of the system R, in the subcritical regime, as function of the number 
of particles n for several values of probability D for the pair-creation conservative contact process. 
The upper straight line has slope 1.33 and the lower one has slope 1. 
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FIG. 10: Log- log plot of size of the system R, in the subcritical regime, as function of the number of 
particles n for several values of probability D for the triplet-creation conservative contact process. 
The upper straight line has slope 1.33 and the lower one has slope 1. 

obtained by Dickman and Tome [s^ show that the transition is continuous for D < 0.95 
which is consistent with our results. However, they argue that the transition should remain 
continuous for any finite diffusion. 

We argue that the formation of a compact cluster {dp = 1) is a signature of a first order 
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transition. First of all, the compact cluster has a nonzero density because p = n/R does not 
vanish in the limit n ^ oo and should therefore be identified with the active phase. Since 
the lattice is infinite the active phase is in coexistence with the frozen state (no particles). 
This behavior is very different from that corresponding to a continuous transition. In this 
case, the fractal dimension is less than one which cannot be identified with the active phase 
since the density p = ~ ^-(i-dF)/c(F _^ q ^^g^ n — > oo. 

When the cluster is a compact one, the ratio n/R gives, in the limit n —>■ oo, the density 
Po of the active phase in coexistence with the frozen phase. We have determined the values 
of Po for several values of D above the tricritical point. An extrapolation for D = 1 gives 
Po = 0.665(1) ^2/3 for the pair-creation model and po = 0.835(2) ^ 5/6 for the triplet- 
creation model. The values of ao at the first order transition can be obtained by substituting 
Po into the exact results ([7j) and ([8]). Using the numerical values, we get ao = 0.222(3) ~ 2/9 
and ao = 0.115(1) ~ 25/216 for the pair-creation model, and the triplet-creation model, 
respectively. As stated before, these values are distinct from the spinodal values coming 
from the exact solutions and ([H]). We remark, on the other hand, that the value of po 
that we have obtained for the triplet-creation model agrees with the value p = 0.84 obtained 
by Dickman and Tome [[8;] for the active coexistence phase at D = 0.95. 

V. CONCLUSION 

The effect of diffusion in nonequilibrium systems has been studied here for the case of three 
conservative contact processes. For the single-creation contact process, the diffusion does 
not change the nature of the phase transition, being continuous for any diffusion rate. This is 
expected since the usual contact process has already an intrinsic diffusion. Indeed, consider 
the following sequence of transitions 010 —>■ Oil — > 001 starting with an isolated particle. 
The net result is a hopping of the isolated particle to a neighboring site, or effectively a 
diffusion. The sequence shown is a possible sequence of states for the single-creation which 
is carried out by a creation followed by a annihilation. This sequence, on the other hand, is 
not possible for the other two models. 

For the pair-creation and and triplet-creation models the transition is continuous for low 
diffusion and becomes discontinuous for high enough diffusion. The present approach in 
which the number of particles is conserved is appropriate to study first order transition 
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because it is possible to distinguish this transition from a continuous one by measuring 
the fractal dimension of the fractal cluster occurring at the critical point. If the fractal 
dimension is smaller than the dimension of the lattice the transition is continuous. When 

the cluster becomes compact, and the fractal dimension equals the lattice dimension, the 
transition becomes first order and, in addition, the density of particles turns out to be the 
density of the active phase in coexistence with the frozen phase. 
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